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Abstract — The pragmatic approach to coded continuous-phase 
modulation (CPM) is proposed as a capacity-achieving low- 
complexity alternative to the serially-concatenated CPM (SC- 
CPM) coding scheme. In this paper, we first perform a selection 
of the best spectrally-efficient CPM modulations to be embedded 
into SC-CPM schemes. Then, we consider the pragmatic capacity 
(a.k.a. BICM capacity) of CPM modulations and optimize it 
through a careful design of the mapping between input bits 
and CPM waveforms. The so obtained schemes are cascaded 
with an outer serially-concatenated convolutional code to form 
a pragmatic coded-modulation system. The resulting schemes 
exhibit performance very close to the CPM capacity without 
requiring iterations between the outer decoder and the CPM 
demodulator. As a result, the receiver exhibits reduced complexity 
and increased flexibility due to the separation of the demodulation 
and decoding functions. 



I. Introduction 

Continuous phase modulation (CPM) is a class of band- 
width efficient modulation schemes [1] whose characteristic 
of constant envelope makes them robust with respect to the 
nonlinearities introduced by the analog baseband and radio- 
frequency sections of low-cost transceivers usually found in 
consumer-type communications equipment. Good application 
examples are the second generation GSM cellular system, and 
satellite communication systems. 

A CPM modulator is a finite-state machine delivering to 
the channel a continuous-phase, constant envelope waveform 
that depends on its input symbol and internal state. In [2] a 
CPM modulator has been shown to be decomposable into the 
cascade of a time-invariant convolutional encoder (continuous- 
phase encoder, CPE) operating on a ring of integers, and 
of a time-invariant memoryless modulator (MM). Recently, 
this decomposition has been exploited by inserting an outer 
convolutional encoder whose coded bits enter an interleaver 
and then the CPE, thus forming what is known in the literature 
as a serially-concatenated convolutional encoder (SCCC) [3]. 
In the following, we will call this scheme SC-CPM. 

Iterating between the outer encoder and the CPE through the 
interleaver yields rather good performance [4], [5], [6], which 
should be compared with the capacity of the CPM scheme. 
Recently, the authors of [7] proposed a simulation-based 
method for the computation of the capacity of channels with 
memory, which can be applied to the evaluation of the capacity 
of the CPM schemes. A similar method has been proposed in 



[8], where it has been applied to a generalized form of CPM 
modulation achieving improved spectral efficiency. 

In communication systems where the channel conditions 
can vary significantly with time, an efficient radio resource 
management requires the availability at the physical layer 
of adaptive coding-modulation, capable of varying its char- 
acteristics of bandwidth and energy efficiency following the 
channel rate of variation. This requirement has originated an 
active research on bit-interleaved (also known as pragmatic) 
coded modulation [9], [10], [11]. Pragmatic coded modulation 
consists in cascading a highly performing versatile binary 
encoder (typically, a punctured turbo or low-density parity- 
check code) with several modulation schemes with increas- 
ingly large signal alphabets. The versatile encoder is capable 
of varying both its rate and codeword length in a wide 
range in order to achieve increasing spectral efficiencies. An 
example of the obtainable results has been published in [12], 
which demonstrates a scheme based on SCCC and linear two- 
dimensional modulations yielding spectral efficiencies in a 
very wide range lying around 1 dB from the Shannon capacity 
limits. This paper deals with the extension of the pragmatic 
approach to CPM modulation (called P-CPM in the following). 
This approach does not require iterations between the outer 
encoder and the CPE, since the CPM is treated exactly as 
a linear modulation in a bit-interleaved turbo-trellis coded 
modulation. A nice consequence is that the overall CPE state 
complexity is not enhanced by the number of iterations, thus 
permitting to increase the bandwidth efficiency through the 
use of a larger number of CPE states. 

In [13] a first attempt to design pragmatic schemes employ- 
ing CPM modulation has been presented. The authors showed 
that the pragmatic capacity of CPM schemes heavily depends 
on the mapping between information bits and CPM signals, 
and presented a simulation-based mapping optimization algo- 
rithm with some examples. Further improvements have been 
presented in [14], where a procedure for optimizing such map- 
ping has been proposed. The optimized CPM modulators have 
been embedded into a P-CPM scheme and its performance has 
been assessed through simulation. For the considered CPM 
schemes, performance improvements of more than 2 dB in 
pragmatic capacity have been observed. 

Previuos literature on the subject includes [15], where 
a scheme consisting of a parallel concatenated turbo code 
and continuous-phase modulation has been investigated and 



DRAFT: June 13, 2008 



2 



a modified encoder has been proposed. In [16], CPM modu- 
lations have been studied in the context of multiple antenna 
systems with layered space-time coding. Reduced-complexity 
receivers have been proposed and differential encoding has 
been introduced in order to obtain an increased coding gain. 

This paper presents a systematic and comprehensive ap- 
proach to the problem of designing capacity approaching SC- 
CPM and P-CPM, through the evaluation of the CPM and 
P-CPM capacities, the search for optimal (in the sense of of- 
fering the best trade-off between bandwidth efficiency, energy 
efficiency and complexity) CPM modulations, their embedding 
into SC-CPM and P-CPM schemes through the optimization 
of the mapping between coded bits and CPM waveforms, and 
a thorough comparison of their error probability performance 
with respect to capacities. Examples refer to rectangular and 
raised cosine frequency waveforms, and to three spectral 
efficiencies deemed important for the applications. 

This paper is organized as follows: in Sec. [TT] a general 
description of the considered coding and modulation systems 
is given. Sec. Ill presents the procedure used for the CPM 
capacity computation through simulation. In Sec. [IV] the 
optimization procedure used to select the best CPM parameters 
to be used in SC-CPM schemes is defined. Sec. [VJpresents the 
procedure used for the optimization of P-CPM schemes: first, 
the optimal mapping is derived, the corresponding optimized 
CPE scheme is given and then the CPM schemes with best 
pragmatic capacity are selected. Sec. VI shows how the 



optimized CPM schemes have been embedded into SC-CPM 



and P-CPM coded systems. Finally, Sec. VII shows the results 
obtained using the selected CPM modulations in SC-CPM and 
P-CPM schemes. 

II. System description 

A CPM modulator is a device with memory that generates 
continuous-phase, constant-envelope modulated waveforms 



x{t) 



2E S 
T 



j </>(«) 



whose phase 



ip(t) = 2ttIi a n q(t — nT) 



(II. 1) 



(H.2) 



depends on the input information symbols a n € 
{±1,±3,...,±(M - 1)}, where M = 2 m is the size 
of the input alphabet. Here, T is the symbol interval, E s is 
the energy per symbol, h = Q/P is the modulation index (Q 
and P are relatively prime integers), and q(t) is the phase 
pulse, a continuous function with the following properties 



<?(*) 



t < 

1 t>LT 



The phase pulse is usually defined as the integral of a 

frequency pulse s(t) 



A CPM scheme is then defined by specifying its parameters 
M, h, L and the frequency pulse s(t). 

In this paper, we will consider rectangular (REC) and raised- 
cosine (RC) frequency pulses. The REC pulse is defined as 



srecO) = o^f K*) 
and the RC pulse is defined as 



i(t - LT)} 



SRc{t) 
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[u(t) - u(t - LT)} 



2nt 



where u(t) is the unitary step function. 

According to the well known Rimoldi decomposition [2], 
the modulator can be represented as the cascade of a con- 
tinuous phase encoder (CPE) and a memoryless modulator 
(MM) as in Fig. [T] The CPE, in general, is a time-invariant 
convolutional encoder operating on a ring of integers. 




Fig. 1. Rimoldi decomposition of the CPM encoder. The block labelled Ep 
is a modulo-P adder. 



The serial structure of the Rimoldi decomposition has been 
exploited by adding an outer convolutional encoder connected 
to the CPE trough an interleaver, so as to form a serially con- 
catenated convolutional encoder (SCCC) with interleaver [17] 
(see Fig. |2]i. The iterative receiver performs decoding iterations 
between the inner SISO decoder, which operates on the CPE 
trellis, and the outer SISO decoder, which operates on the outer 
convolutional code trellis. This way, very good performance 
can be achieved, to be compared with the capacity of the 
CPM scheme. In the following, we will denote this scheme 
as serially concatenated CPM (SC-CPM) [4]. 

In this paper, we also consider a pragmatic approach, which 
results in a coded modulation scheme equivalent to a bit- 
interleaved coded modulation (BICM, see Fig.[3]l. It consists in 
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Fig. 2. Block diagram of a serially concatenated CPM co-decoder. The block 
labeled Ai maps blocks of m bits to M-ary CPM symbols. 
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Fig. 3. Block diagram of an encoded CPM mo-demodulator used with a 
pragmatic approach. 



cascading a powerful variable-rate binary encoder (a serially- 
concatenated convolutional code (SCCC) with interleaver in 
this paper) with the CPM modulator. The CPM modulator 
input is connected to the output of the channel encoder 
that computes the binary sequence b = (bi,i G Z) from 
the binary information sequence d. Then,the coded binary 
sequence is mapped to the CPM modulator that generates a 
corresponding CPM signal sequence x(t; b). We will call this 
scheme pragmatic CPM (P-CPM). 

The channel is an additive white Gaussian noise (AWGN) 
channel, whose output signal is y(t; b) = x(t; h)+n(t), where 
n(t) is a zero-mean white Gaussian process with two-sided 
power spectral density No/2. At the receiver, the CPM soft 
demodulator and demapper provide the sequence A = (A, , i £ 
Z) of log-likelihood ratios on the information bit sequence to 
the outer decoder 



A, = log 



P(bj = l\y(t;b)) 
P(bi = 0\y(t;b)) 



Finally, A is used by the outer iterative decoder to compute 
the information sequence estimate d. 

The main difference between P-CPM and SC-CPM is that 
no iterations are performed between the inner CPM demodula- 
tor and the outer turbo decoder. Thus, it keeps the nice features 
of all pragmatic approaches, which merge independent binary 
codes with higher order modulations without requiring joint 
optimization. 

Our goal is to design coded modulation schemes based on 
CPM capable of transmitting information at rates close to 
the CPM channel capacity^] Such goal will be pursued by 
analyzing the P-CPM and SC-CPM schemes and choosing the 
best trade-off between the following characteristics: 

• Spectral efficiency. 

• Energy efficiency. 

• Decoding complexity. 

Spectral efficiency and energy efficiency are two well- 
known concepts. The decoding complexity is defined here as 
the overall number of trellis edges per information bit visited 
by the decoding algorithm. This definition depends on the type 

1 Hereafter, with CPM channel capacity we mean constrained capacity with 
uniform input distribution. 



of coded modulation being used (i.e., P-CPM or SC-CPM). For 
the P-CPM scheme, we have 



J^P-CPM = J^SCCC + !VcPM 



(H.3) 



where 



}cpm — MN S CPE 



No 
mK n 



and N Q and K a are, respectively, the outer code word length 
and information word length, and N s ,cpe = PM L ~ l is the 
number of states of the CPE. The complexity of the SCCC 
binary decoder is J^sccc an d will be defined later according 
to the chosen code (see dVILTT i), 

For the SC-CPM scheme, we adopt the following definition, 
which holds for binary convolutional codes with rate ^f- 
obtained by puncturing a rate 1/2 mother code: 



3^sc-CPM — Nn (2N so + J^cpm) 



(H.4) 



where N t is the number of decoding iterations between the 
outer SISO and the inner CPE SISO, and N so is the number 
of states of the outer convolutional encoder. 

In the following, we will optimize the CPM schemes to be 
embedded into both SC-CPM and P-CPM starting from their 
capacity evaluation. 

III. CPM CAPACITY COMPUTATION 

With reference to Figures [2] and [3] CPM signals are infinite- 
length waveforms. Consider then a finite observation window 
[—NT, NT] of 2^+1 symbol intervals, and define the channel 
mutual information over it: 



Z(X;Y)= J B x , y hog 2 



Pv|x(y|x) 
pr(y) 

where X is a vector of samples of x(t) that form a sufficient 
statistic of x(t) in the interval [—NT, NT] and Y is the 
correspondent vector of channel outputs. The set of all possible 
values for X will be denoted X throughout the paper. The 
CPM capacity can be defined through the limit 



a 



CPM 



n^L (2N+l)T 



7(X;Y) [bits/s] (III.l) 



Using the definition of mutual information I and of the 
log-likelihood ratio (LLR) 

A(u,y) = l0g 2 (py|x(y|u)) -log 2 (p Y |x(y|Uref)) 

as well as the max* operator [17] 

max* (a, 6) = log 2 (2 a + 2 b ) 
we can transform ( |III.l| i into 
1 



a 



CPM 



= to— lim 



i^o (2N + l)f E *>» \^ a / A(u ' y) - A(x ' y) 

where to is the number of bits per CPM symbol. 

It has been observed in [7] that the first term inside the 
average is a by-product of the SISO algorithm and corresponds 
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to the max* of the forward path metrics a at step N in the 
CPE decoder 

max*A(u,y) =max*a w (s), 

u£X s 

with the following initialization 

0, a = 0; 



a_jv(s) = 



— oo, s =/= 0. 



where s runs in the set of trellis states of the CPE. Moreover, 
since the channel is memoryless, the second term is obtained 
by summing the LLRs of the transmitted waveforms: 



A ( X 'V) = E X ( X iiVi)- 



i=-N 



Finally, invoking the ergodic properties of the system, the 
ensemble average can be removed leading to: 



Ccpm = Tn — lim 



max € X*X(u, y) — A(x, y) 



n^oc 2N + 1 V 

Thus, the CPM capacity can be estimated through a Monte- 
Carlo simulation of the three internal blocks of Fig. [3] i.e., the 
CPM modulator, the channel, and the CPM soft demodulator, 
followed by a time average. 

IV. Optimization procedure 

We define the complexity of a CPM scheme as the number 
of edges per CPM input bij^] y = P2 mL /m, where the 
parameters have been defined in Section II. Our aim is to 
maximize the CPM capacity measured in bits/s/Hertz versus 
the signal-to-noise ratio for a given complexity. The optimiza- 
tion algorithm modifies the CPM parameters (to, P, and L) 
yielding the given complexity, and computes the CPM capacity 
expressed in bits per CPM waveform for a given signal-to- 



noise ratio according to the algorithm described in Sec. Ill 
Then, based on a bandwidth definition of the CPM scheme, it 
evaluates the rate and chooses the best scheme for all signal- 
to-noise ratios of interest. In the following, we describe the 
optimization algorithm step by step: 

1) Define the bandwidth B of the CPM system that contains 
a given percentage of the total signal power. 

2) Define the CPM symbol rate as R s = l/T. 

3) The symbol signal-to-noise ratio E s /Nq is given by 

E. = J_Pt 
N a R s N 

where Pt is the transmitted power. 

4) Compute the capacity (in bits/symbol) of each CPM 
scheme using the method described in Section [TTT] as a 
function of the symbol SNR Ccpm = Ccpm(E s /N ). 

5) Evaluate the CPM capacity as 



C ReC( 



CPM 



1 P T 
RsNo 



[bits/s/Hz] 



(IV. 1) 



2 Please note that the definition of y differs from 3^cpm m that the former 
is the number of CPE trellis edges per CPM input bit, while the latter is the 
number of CPE trellis edges per bit at the channel encoder input. 



6) All CPM schemes with a given complexity y are com- 
pared with respect to C for each value of j^P- and the 
best is chosen. 

7) Finally, to obtain the normalized plot that uses the bit 
SNR Eb/No on the abscissa we use the relationship 

Eh _ 1 E s _ R s E s _ 1 Pt 

~ ~ ~C iVo " 



(IV.2) 



N Ccpm N C N C N 
In Fig. [4] and Fig.[5]we show the capacity versus E^/Nq for 
the best CPM schemes with REC and RC frequency pulses 
and complexity ranging from 8 to 512. The bandwidth is 
defined as the one including 99% of the total signal power. 
For comparison purposes, in the figure we have also plotted 
the capacity of QPSK and 8PSK modulations with a square 
root raised cosine shaping filter with roll-off 0.25, and the 
unconstrained capacity of the AWGN channel. 

In Table [I] and Table |5| we list for each SNR (E S /N Q ) 
and complexity from 8 to 512 the CPM schemes with REC 
and RC frequency pulses achieving the highest capacity for a 
given signal-to-noise ratio. Each scheme is characterized by 



the parameters E{,/Nq as defined in ( IV.2 1, the capacity C as 
defined in ( IV.l} , the CPM parameters to, L, P, and finally 
the symbol rate R s . 




e t //v„ [dB] 

Fig. 4. Capacity versus E^/Nq of the best CPM schemes with rectangular 
frequency pulse and variable complexity. The bandwidth is defined at 99 % 
of the total power. 



V. Pragmatic capacity of a TCM scheme 
With reference to Fig. [3] let the binary sequence b be formed 
by vectors B n = (B n ,i, • • • , B njll ) of the to bits entering 
the CPM modulator at the n-th trellis step. As previously, we 
consider a length-(2A^ + 1) input sequence of B„. 

The scheme of Fig. [3] can be characterized by a pragmatic 
capacity Cp_cpm> defined by: 



Cp_cpm - — 



N m 

T E E 

n=-N i=l 



Y), (VI) 



which depends on the mapping between input binary se- 
quences and CPM waveforms. 

In order to optimize the system performance, it is crucial to 
search for the optimal mapping, i.e., the one that maximizes 
the pragmatic capacity of a given CPM scheme. To this 
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TABLE I 



Table of the CPM schemes with best C and rectangular frequency pulse at 99% bandwidth. The complexity range is 8 to 5 1 2. The 

HIGHLIGHTED ENTRIES CORRESPOND TO THE SELECTED CPM SCHEMES FOR TARGET CAPACITIES OF 1, 1.5 AND 2 BITS/S/HZ. 
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TABLE II 



Table of the CPM schemes with best C and raised-cosine frequency pulse at 99% bandwidth. The complexity range is 8 to 512. 
The highlighted entries correspond to the selected CPM schemes for target capacities OF 1, 1.5 AND 2 bits/s/Hz. 



purpose, we will derive in the following a lower bound to the 
pragmatic capacity, which makes explicit the dependence on 
the mapping. Before delving deeper into the analysis, though, 
we show in Fig. [6] a comparison between the optimized P- 
CPM capacity, obtained through the algorithm derived in the 
following, and that obtained by a straightforward application 
of the mapping induced by the Rimoldi decomposition. The 
curves show that the pragmatic capacity of the optimized CPE 
is significantly improved with respect to the non optimized 
CPE. For the non optimized CPE, the gap between the CPM 
capacity and the pragmatic capacity is roughly 1.5 dB, while 
the pragmatic capacity of the optimized CPE is very close to 
the CPM capacity. The obtained gain is roughly 1.4 dB. 



Since 

J(B Bii ;Y) = l-H(B n>i \ Y), (V.2) 

the optimal mapping is the one that yields a minimum of 
Eft£<#0Bn,<|Y). Define X n .,{b), 6-0,1, as the set of 
CPM waveforms whose mapping satisfies B n i = b. Obvi- 
ously, X = X n ,i(0) U Xn,i{\). We have: 

H(B n ,i\Y) = J H{B n>i \Y = y)p Y (y)rfy. (V.3) 
Now, H(B n: i\Y — y) is given explicitly by: 

g(s „„iY =y)=fl f £ "^'°"'7; y i u) y (v.4) 
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E t /W„ [dB] 

Fig. 5. Capacity versus E^/Nq of the best CPM schemes with raised-cosine 
frequency pulse and variable complexity. The bandwidth is defined at 99 % 
of the total power. 




CPM capacity 
Pragmatic capacity (opt.) 
Pragmatic capacity (non opt.) 



2 3 4 

E„IN„ [dB] 



Fig. 6. CPM capacity and pragmatic capacities of the binary REC scheme 
with h = 1/2 and L = 3. 



By Jensen's inequality and the concavity of the entropy 
function, we obtain the upper bound: 



H(B n<i \Y) < i E W 

6=o,i y 



H(B nti \Y) <\Y. H \ 

6=0,1 V 



2_ Eue^.(o) Py|x(x]u ) 

(V.7) 

Owing to the ergodic property of the system, it can be shown 
that EugA" Py|x( x |u) is equal to a constant for almost every 
transmitted x. Thus: 

' 2 E x g.y„, t (fc) Sue* n , ( (0) Py|x(x|u) 

Ex£A'EueA'PY|x(x|u) 

(V.8) 

Finally, notice that, since Py|x( x I u ) = Py|x( u I x ) f° r the 
AWGN channel, the two entropies in the RHS of ( |V.8| > have 
the same arguments, and thus are equal: 

' 2 (o) Eu6^, 4 (0) Py|x(x|u) 1 



H(B n ,i\Y)<H 



£xe*£ue*PY|x(x|u) 



(V.9) 



or, in a more explicit form: 



H(B nti \Y)<- 



2 E, 



xeAT„,i(0) EugAr„,i(0) ^Y| 

ExeA-Eue;tPY|x(x|u) 
b 2 Exe^o) E u gA-„.,(o) Py|x(x|u) 
ExeArEueA-PY|x(x|u) 

2 ExgA' n , i (0)EugA', 1 , i (l)PY|x(x|u) 

ExeA-EuGA-PYIxlxlu) 
b 2 ExgA-„.,(0) EugA-„.,(l) PY| X (x|u) 

Exe*Eu e *PY|x(x|u) 
The above equation, substituted in dVZk and, then, back in 



Pyix(x|u) 



(V.l i, gives the lower bound on the pragmatic capacity that is 
the starting point for our optimization. 



where the binary entropy function H(p) = — plogp — (1 — 
p) log(l — p) has been used in the RHS of the above equation. 

To find a viable path to mapping optimization, we approxi- 
mate the distribution of Y with its value for infinite SNR, by 



making in (V.3i the following substitution: 



1 

\x\ 



(V.5) 



where \X\ is the size of X and S(-) is Dirac's delta. The above 
approximation makes sense because the optimal mapping for 
large SNR should reasonably be the same as for infinite SNR. 



Substituting (V.5i into (V.3i, we obtain: 



H(B n ,\Y)J- E 



2 ^ \X\ 

6=0,1 1 1 



E » 

xe^„,j(6) 



EuE^„,,(0) PY|x(x|u)' 

E ue ^PY|x(x|u) 



(V.6) 



where the sum w.r.t x over X has been split into separate sums 
over X n}i (0) and X nii (l). 



A. Mapping optimization 

The problem of mapping optimization has been reduced 



to the minimization of the upper bound in (V.9i (or (fv])). 
A binary entropy function is minimized if its argument dis- 
tribution is made as unbalanced as possible. For large SNR, 
Exe* (o) E ue *„, ( (o) PY|x(x, u) will be the dominant term 
in ( jV.9[ ), because it contains the term with x = u; thus, 
it should be made as close to 1 as possible. Conversely, 
the complementary term E x eAf n>j (o) Eue* n>i (i) Py|x(x, u), 
should be made as close to as possible, i.e., it should be 
minimized. We have come up with the following design rule 
for the optimal mapping: 



<*n,i(0) = arg 



mm 

X' C X 
\X'\ = \X\/2 x£X ' ™$ x ' 



E E PY|x(x,u). 



(V.10) 



For practical reasons, the mapping is generated by an m- 
bit labelling of the trellis edges. Let 7^(0) (7^(1)) be the 
subset of trellis edges at time n whose label has a (1) in 
the z-th position. Then, X n i(0) is constituted by all CPM 
waveforms whose trellis paths pass through an edge belonging 



to 7^,(0), and the design rule in (V.lOi can be restated in 
terms of 7^(0). 
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We impose the condition that the trellis labelling must be 
right-resolving, i.e., edges leaving the same trellis state have 
different binary labels. This condition does not force a trellis 
labelling that varies with n. Moreover, the design criterion in 



(V.lOi is also independent of n. From these facts, we deduce 
that the optimal mapping is generated by a time-invariant 
trellis labelling. Thus T n ^(0) = %(0), for every n. Instead, 
the condition of having a right-resolving trellis implies that 
7i(0), . . . , %n(0) cannot be chosen independently, for m > 1. 

In practice, we approximate the design rule in ( V. 10| > 
with the following suboptimal rule, which derives from the 
assumption that, for large SNR, there is one dominant term in 
the sum of ( |V.10) : 



*n,i(0) 



are mm 

X'CX 



max max Pyix( x > u ) 
xex'ugx' 1 



(V.ll) 



Based on ( | V. 1 1| >, our mapping optimization procedure consists 
then of the following steps: 

* For every pair of trellis edges, we apply the BCJR 
algorithm to compute max u max x pY|x( x , u ) over all 
pairs of paths passing through the given pair of edges 
at time zero. To do this in general, we have to extend the 
trellis both at the left (negative time instants) and at the 
right (positive time instants) of the zeroth section. After 
a few trellis steps, the path metrics reach a steady-state 
value, so there is no need to proceed further. 

• The previous step yields a metric for each pair of edges. 
Since we want to minimize the expression between braces 

we partition the edges into M = 2" 



in (V.ll 



equal-size 

clusters, corresponding to the M different ?7i-bit labels, 
in such a way that pairs of edges with the highest metric 
are all clustered together. In doing this, we force pairs of 
edges with the same starting state into different clusters, 
to allow for a right-resolving labelling. 
• We map the M = 2 m clusters to m-tuples of bits accord- 
ing to a Gray mapping. More precisely, we define the 
distance between two clusters as the minimum distance 
between CPM waveforms passing at time zero through 
edges belonging to the two clusters. Cluster pairs with 
the smallest distance will be associated to binary labels 
with Hamming distance 1. 
The highest metric corresponds to edge pairs that belong to 
the same trellis section of a minimum-distance error event. It 
may be questionable whether it is possible to cluster all pairs 
of edges with that property, apart from those with the same 
starting state, and whether this gives a unique clustering. In 
the next subsection, we will completely answer these questions 
for M — 2, and give the analytical expression of the optimal 
mapping, provided that some necessary conditions are met. 

B. The optimal mapping for M = 2 

As described in the previous subsection, the algorithm for 
mapping optimization puts into the same cluster all edge pairs 
that belong to the same trellis section of a minimum-distance 
error event. For labelling purposes, we would like to obtain in 
this way exactly M clusters, but this happens only at certain 
conditions. In this subsection, we derive such conditions for 
binary CPM schemes. 



Let us consider a CPM scheme with M — 2, impulse 
length L and modulation index h = Q/P. Trellis edges will 
be denoted hereafter through the couple (a,f3), where a is a 
length-!/ binary vector, including the correlative state and the 
input symbol, and j3 G {0, . . . , P — 1} represents the phase 
state. Let xi(t) and X2(t) be two CPM waveforms, corre- 
sponding to input symbol sequences ai and a2, respectively. 
It is well known that the Euclidean distance between these 
two CPM waveforms only depends on the difference sequence 
b = ai &2, with elements belonging to {—1, 0, 1}. 

Now, consider a given difference sequence b and define a 
graph G(h) = (V,£(b)), where: 

• the vertex set V is in one-to-one correspondence to the 
set of trellis edges for the considered CPM scheme. 

• the edge set £(b) is constructed in the following way: 
two vertices are connected by an edge if and only if the 
corresponding trellis edges have different starting states 
and belong to the same trellis section of an error event 
generated by b. 

Let Co, . . . ,Cjv(b)-i be the Nih) connected components 
of the graph £7(b). The following theorem gives the important 
properties of C , . • . ,C^( b )_i. 

Theorem 5.1: Consider a CPM scheme with M = 2, 
impulse length L and modulation index h = Q/P- Given a 
difference sequence b, with length A(b), the following facts 
about the graph GO 3 ) defined above hold: 

1) If A(b) = 2, then Af(b) = P, otherwise W(b) = 1. 

2) If JV(h) = P, the connected component Ci, i — 
0, . . . , P — 1 is the subgraph induced by the subset of 
vertice^O 

Vi = {(a, (« - Qw H (a)) P ) : a G {0, 1} L } , (V.12) 

independently of b. All C/s have size 2 L . 
Proof: See Appendix [I] ■ 

Suppose now that b gives the minimum Euclidean distance 
of the CPM scheme. We want to partition the edges into two 
clusters, namely, T(0) and T(l) (edges with a label and 
edges with a label 1, respectively). To meet all constraints 
given by the graph <?(b), every subset Vi should be entirely 
contained into one of the two clusters. This is clearly not 
possible if A(b) ^ 2, because in that case Vi = V. 

If A(b) = 2, instead, we can place exactly P/2 subsets 
into T(0) and the other P/2 into T(l). (Notice that the sizes 
of T(0) and of T{1) both must be equal to P2 L ~ 1 .) This can 
be done only if P is even, otherwise, one of the Vi's must be 
split in two. Thus, we have proved the following proposition: 

Proposition 5.1: Consider a CPM scheme with M = 2, 
impulse length L and modulation index h = Q/P. We can 
fully perform the algorithm of mapping optimization described 
in the previous section if and only if: 

• Every difference sequence b giving the minimum Eu- 
clidean distance of the scheme has length two, and 

• P is even. 

3 The notation lp means "I modulo P". 
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If the conditions of Prop. |5.l| are satisfied, then, the optimal 
mapping is uniquely determined by the condition of right- 
resolving labelling. In fact, the two trellis edges leaving a given 
trellis state will belong to Vi and V( i+ g) p , for some i. Thus, 
for the labelling to be right-resolving, it must be: 



and 



T(0) = V U V [2Q)p U • • • U V ((P _ 2 )Q )p (V.13) 



T(l) = V Qp U V (3Q)p U ■ • • U V ((P _ 1)Q)p , (V.14) 



or vice versa. In this way, all CPM waveform pairs at a 
minimum distance between each other will be associated to 
binary input sequence pairs with Hamming distance equal to 
1. 

If the conditions of Prop. |5.1| are not satisfied, then the 
mapping optimization algorithm of the previous section cannot 
be fully performed. It means that some CPM waveform pairs 
with minimum distance will be associated to binary input 
sequence pairs with Hamming distance larger than 1. A 
reasonable approach in this case is to try to minimize the 
number of such CPM waveform pairs. 

C. The optimal mapping for M > 2 

In non-binary CPM schemes, things become more involved. 
In general, we have not found necessary and sufficient condi- 
tions to fully perform the algorithm of mapping optimization. 



Here, we give a generalization of Theorem 5.1 which how 




ever deals only with a subset of possible difference sequences. 

Theorem 5.2: Consider a CPM scheme with M > 2, 
impulse length L and modulation index h = Q/P. Let b = 
(&i, . . . , &A(b)) be a difference sequence with length A(b) 
and elements in {0, ±1, . . . , ±(M — 1)}. Let &A(t>) = ±L 
The following facts about the graph 5(b) defined above hold: 

1) If A(b) = 2, then Af(b) = P, otherwise W(b) = 1. 

2) If Af(h) = P, the connected component C p , p — 
0, . . . , P — 1 is the subgraph induced by the subset of 
vertices 



(V.15) 

a=( ai ,...,a L )£{0,...,M-l} L }, 

independently of b. All C p 's have size M L . 
Proof: See Appendix |lT| ■ 

Following the same arguments of the previous subsection, 
we thus obtain the following proposition, which gives only 
sufficient (not necessary) conditions to fully perform the 
mapping optimization algorithm. 

Proposition 5.2: Consider a CPM scheme with M > 2, 
impulse length L and modulation index h = Q/P. We can 
fully perform the algorithm of mapping optimization described 
in the previous section if: 

• Every difference sequence b giving the minimum Eu- 
clidean distance of the scheme has length two and 6 2 = 
±1, and 

• P is a multiple of M. 



If the conditions of Prop. |5.2| are satisfied, then, by con- 
straining the resulting labelling to be right-resolving, we obtain 
the following M clusters: 

^(*) = V(iQ)p U V((M+i)Q) P U • • • U V((P-M+i)Q) P , (V.16) 

for i = 0, . . . , M - 1. 

Finally, the clusters T(i)'s are mapped to the M binary 
labels according to Gray mapping. 

In the cases that are not covered by the hypotheses of 



Theorem 5.2 the algorithm may or may not be fully performed 
depending on the scheme parameters and on the difference 
sequences b giving the minimum Euclidean distance. If the 
algorithm can be fully performed, in general, the M clusters 



may not be as in (V. 16 1 



D. The optimized CPE 



Starting from (V.15 1, where we define 




we derive the input symbol a\ to the CPE as the solution to 
the following equation 



(Qai)p = [p - a L+1 



1=2 / 



Since P and Q are relatively prime, such solution is unique. 
Moreover, from (V.15 i we have a\ <E {0, . . . ,M — 1} and, 
since also M and Q are relatively prime, (QojJm entails a 
unique solution for As a result, we can write 



(Qai 



M 




(V.17) 



Clustering as in (V. 16 1 can be performed observing that 
the vertex subset indices in each cluster can be written as 



p = ((kAl + i)Q)p. Inserting this in ( V. 1 7 1 and recalling that 
P is a multiple of M we obtain 



{Qa\)M = iQ + kMQ - a L+1 - Q ^ a ; 

V 1=2 / 

= [iQ - a L+1 - Qy^a; 

V 1=2 J , 



where i is the cluster label. Applying the distributive property 
of multiplication over addition in the ring Zm we obtain 



a i 



(a-L+iQ V 



M 




(V.18) 



Finally, the binary label B is obtained as B = G(i), where G 
is the Gray function defined in [18, ch.20]. Fig. [7] shows the 
scheme of the optimized CPM encoder. 

Adopting the optimized CPE, i.e., starting from the CPE 
structure of Fig. [7] a procedure similar to the one described 
in Sec. IV has been performed in order to select the CPM 
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1.28 


2 1 


4 


2 


0.96 


3.92 


1.28 


2 1 


4 


2 


0.96 


3.92 


1.28 


2 1 


4 


2 


0.96 


3.92 


1.28 


2 1 


4 


2 


0.96 


3.92 


1.28 


2 1 


4 


2 


0.96 


5.5 


4.64 


1.22 


2 1 


4 


1 


0.65 


4.43 


1.28 


2 1 


5 


1 


0.8 


4.25 


1.34 


2 1 


4 


2 


0.96 


4.25 


1.34 


2 1 


4 


2 


0.96 


4.25 


1.34 


2 1 


4 


2 


0.96 


4.25 


1.34 


2 1 


4 


2 


0.96 


4.25 


1.34 


2 1 


4 


2 


0.96 


6 


5.07 


1.24 


2 1 


4 


1 


0.65 


4.77 


1.33 


2 1 


5 


1 


0.8 


4.58 


1.39 


2 1 


4 


2 


0.96 


4.58 


1.39 


2 1 


4 


2 


0.96 


4.58 


1.39 


2 1 


4 


2 


0.96 


4.58 


1.39 


2 1 


4 


2 


0.96 


4.58 


1.39 


2 1 


4 


2 


0.96 


6.5 


5.52 


1.25 


2 1 


4 


1 


0.65 


5.13 


1.37 


2 1 


5 


1 


0.8 


4.93 


1.44 


2 1 


4 


2 


0.96 


4.93 


1.44 


2 1 


4 


2 


0.96 


4.93 


1.44 


2 1 


4 


2 


0.96 


4.93 


1.44 


2 1 


4 


2 


0.96 


4.93 


1.44 


2 1 


4 


2 


0.96 


7 


5.97 


1.27 


2 1 


4 


1 


0.65 


5.51 


1.41 


2 1 


5 


1 


0.8 


5.29 


1.48 


2 1 


4 


2 


0.96 


5.29 


1.48 


2 1 


4 


2 


0.96 


5.29 


1.48 


2 1 


4 


2 


0.96 


5.29 


1.48 


2 1 


4 


2 


0.96 


5.29 


1.48 


2 1 


4 


2 


0.96 


7.5 


6.44 


1.28 


2 1 


4 


1 


0.65 


5.90 


1.45 


2 1 


5 


1 


0.8 


5.64 


1.54 


1 1 


4 


3 


2.29 


5.64 


1.54 


1 1 


4 


3 


2.29 


5.64 


1.54 


1 1 


4 


3 


2.29 


5.64 


1.54 


1 1 


4 


3 


2.29 


5.64 


1.54 


1 1 


4 


3 


2.29 


8 


6.92 


1.28 


2 1 


4 


1 


0.65 


6.31 


1.48 


2 1 


5 


1 


0.8 


5.97 


1.60 


1 1 


4 


3 


2.29 


5.95 


1.60 


1 1 


4 


4 


2.72 


5.95 


1.60 


1 1 


4 


4 


2.72 


5.95 


1.60 


1 1 


4 


4 


2.72 


5.95 


1.60 


1 1 


4 


4 


2.72 


8.5 


7.40 


1.29 


2 1 


4 


1 


0.65 


6.72 


1.51 


2 1 


6 


1 


0.88 


6.31 


1.66 


1 1 


4 


3 


2.29 


6.27 


1.67 


1 1 


4 


4 


2.72 


6.27 


1.67 


1 1 


4 


4 


2.72 


6.27 


1.67 


1 1 


4 


4 


2.72 


6.27 


1.67 


1 1 


4 


4 


2.72 


9 


7.89 


1.29 


2 1 


4 


1 


0.65 


7.07 


1.56 


2 1 


7 


1 


0.96 


6.66 


1.72 


1 1 


4 


3 


2.29 


6.59 


1.74 


1 1 


4 


4 


2.72 


6.59 


1.74 


1 1 


4 


4 


2.72 


6.59 


1.74 


1 1 


4 


4 


2.72 


6.59 


1.74 


1 1 


4 


4 


2.72 


9.5 


8.38 


1.29 


2 1 


4 


1 


0.65 


7.43 


1.61 


2 1 


7 


1 


0.96 


7.02 


1.77 


1 1 


4 


3 


2.29 


6.92 


1.81 


1 1 


4 


4 


2.72 


6.92 


1.81 


1 1 


4 


4 


2.72 


6.92 


1.81 


1 1 


4 


4 


2.72 


6.92 


1.81 


1 1 


4 


4 


2.72 


10 


8.88 


1.30 


2 1 


4 


1 


0.65 


7.79 


1.66 


2 1 


8 


1 


1.03 


7.39 


1.82 


1 1 


4 


3 


2.29 


7.27 


1.88 


1 1 


4 


4 


2.72 


7.27 


1.88 


1 1 


4 


4 


2.72 


7.27 


1.88 


1 1 


4 


4 


2.72 


7.27 


1.88 


1 1 


4 


4 


2.72 


10.5 


9.38 


1.30 


2 1 


4 


1 


0.65 


8.15 


1.72 


2 1 


8 


1 


1.03 


7.77 


1.88 


1 1 


4 


3 


2.29 


7.63 


1.94 


1 1 


4 


4 


2.72 


7.63 


1.94 


1 1 


4 


4 


2.72 


7.63 


1.94 


1 1 


4 


4 


2.72 


7.63 


1.94 


1 1 


4 


4 


2.72 


11 


9.87 


1.30 


2 1 


4 


1 


0.65 


8.53 


1.77 


2 1 


8 


1 


1.03 


8.16 


1.93 


1 1 


4 


3 


2.29 


7.99 


2.00 


1 1 


4 


4 


2.72 


7.99 


2.00 


1 1 


4 


4 


2.72 


7.99 


2.00 


1 1 


4 


4 


2.72 


7.99 


2.00 


1 1 


4 


4 


2.72 


11.5 


10.37 


1.30 


2 1 


4 


1 


0.65 


8.91 


1.82 


2 1 


8 


1 


1.03 


8.55 


1.97 


1 1 


4 


3 


2.29 


8.36 


2.06 


1 1 


4 


4 


2.72 


8.36 


2.06 


1 1 


4 


4 


2.72 


8.36 


2.06 


1 1 


4 


4 


2.72 


8.36 


2.06 


1 1 


4 


4 


2.72 


12 


10.87 


1.30 


2 1 


4 


1 


0.65 


9.31 


1.86 


2 1 


8 


1 


1.03 


8.95 


2.02 


1 1 


4 


3 


2.29 


8.73 


2.12 


1 1 


4 


4 


2.72 


8.73 


2.12 


1 1 


4 


4 


2.72 


8.73 


2.12 


1 1 


4 


4 


2.72 


8.73 


2.12 


1 1 


4 


4 


2.72 


12.5 


11.37 


1.30 


2 1 


4 


1 


0.65 


9.72 


1.90 


2 1 


8 


1 


1.03 


9.36 


2.06 


1 1 


4 


3 


2.29 


9.12 


2.18 


1 1 


4 


4 


2.72 


9.12 


2.18 


1 1 


4 


4 


2.72 


9.12 


2.18 


1 1 


4 


4 


2.72 


9.12 


2.18 


1 1 


4 


4 


2.72 


13 


11.87 


1.30 


2 1 


4 


1 


0.65 


10.14 


1.93 


2 1 


8 


1 


1.03 


9.78 


2.10 


1 1 


4 


3 


2.29 


9.50 


2.24 


1 1 


4 


4 


2.72 


9.50 


2.24 


1 1 


4 


4 


2.72 


9.50 


2.24 


1 1 


4 


4 


2.72 


9.50 


2.24 


1 1 


4 


4 


2.72 


13.5 


12.37 


1.30 


2 1 


4 


1 


0.65 


10.58 


1.96 


2 1 


8 


1 


1.03 


10.20 


2.14 


1 1 


4 


3 


2.29 


9.89 


2.29 


1 1 


4 


4 


2.72 


9.89 


2.29 


1 1 


4 


4 


2.72 


9.89 


2.29 


1 1 


4 


4 


2.72 


9.89 


2.29 


1 1 


4 


4 


2.72 


14 


12.87 


1.30 


2 1 


4 


1 


0.65 


11.02 


1.99 


2 1 


8 


1 


1.03 


10.64 


2.17 


1 1 


4 


3 


2.29 


10.29 


2.35 


1 1 


4 


4 


2.72 


10.27 


2.36 


1 1 


6 


4 


3.61 


10.27 


2.36 


1 1 


6 


4 


3.61 


10.27 


2.36 


1 1 


6 


4 


3.61 


14.5 


13.37 


1.30 


2 1 


4 


1 


0.65 


11.48 


2.01 


2 1 


8 


1 


1.03 


11.08 


2.20 


1 1 


4 


3 


2.29 


10.69 


2.40 


1 1 


4 


4 


2.72 


10.62 


2.45 


1 1 


6 


4 


3.61 


10.62 


2.45 


1 1 


6 


4 


3.61 


10.62 


2.45 


1 1 


6 


4 


3.61 


15 


13.87 


1.30 


2 1 


4 


1 


0.65 


11.94 


2.02 


2 1 


8 


1 


1.03 


11.53 


2.22 


3 1 


10 


1 


0.77 


11.08 


2.47 


2 1 


8 


2 


1.61 


10.97 


2.53 


1 1 


6 


4 


3.61 


10.97 


2.53 


1 1 


6 


4 


3.61 


10.97 


2.53 


1 1 


6 


4 


3.61 



TABLE III 



Table of the P-CPM schemes with best Cp_cpm and rectangular frequency pulse at 99% bandwidth. The complexity range is 8 to 
512. The highlighted entries correspond to the selected CPM schemes for target capacities of 1, 1.5 and 2 bits/s/Hz. 



Es/NO 


8 


16 


32 


64 


128 


256 


512 


Lb NO 


Cp rr 


P 


L 


Rs 


Lb NO 


Cp 


Tl Q 


P 


L 


Rs 


Lb NO 


Cp 


n Q 


P 


L 


Rs 


Lb \u 


Cp 


n Q 


P 


L 


Rs 


Lb NO 


Cp 


n Q 


P 


L 


Rs 


Lb Lib 


Cp 


n Q 


P 


L 


Rs 


Eb/NO 


Cp 


n Q 


P 


L 


Rs 





1.52 


0.71 


1 2 


2 


0.94 


1.16 


0.77 


1 1 


2 


3 


1.14 


0.96 


0.80 


1 1 


2 


4 


1.31 


0.93 


0.81 


1 1 


2 


5 


1.44 


0.93 


0.81 


1 1 


2 


5 


1 .44 


0.93 


0.81 


1 1 


2 


5 


1.44 


0.93 


0.81 


1 1 


2 


5 


1.44 


0.5 


1.81 


0.74 


1 2 


2 


0.94 


1.42 


0.81 


1 1 


2 


3 


1.14 


1.20 


0.85 


1 1 


2 


4 


1.31 


1.14 


0.86 


1 1 


2 


5 


1.44 


1.14 


0.86 


1 1 


2 


5 


1.44 


1.14 


0.86 


1 1 


2 


5 


1.44 


1.14 


0.86 


1 1 


2 


5 


1.44 


1 


2.12 


0.77 


1 2 


2 


0.94 


1.69 


0.85 


1 1 


2 


3 


1.14 


1.43 


0.91 


2 1 


4 


2 


0.8 


1.38 


0.92 


1 1 


2 


5 


1.44 


1.38 


0.92 


1 1 


2 


5 


1.44 


1.38 


0.92 


1 1 


2 


5 


1.44 


1.38 


0.92 


1 1 


2 


5 


1.44 


1.5 


2.46 


0.80 


1 2 


2 


0.94 


1.99 


0.89 


1 1 


2 


3 


1.14 


1.64 


0.97 


2 1 


4 


2 


0.8 


1.63 


0.97 


1 1 


2 


5 


1.44 


1.63 


0.97 


1 1 


2 


5 


1.44 


1.63 


0.97 


1 1 


2 


5 


1.44 


1.63 


0.97 


1 1 


2 


5 


1.44 


2 


2.81 


0.83 


1 2 


2 


0.94 


2.30 


0.93 


1 1 


2 


3 


1.14 


1.86 


1.03 


2 1 


4 


2 


0.8 


1.86 


1.03 


2 1 


4 


2 




1.86 


1.03 


2 1 


4 


2 


0.8 


1.86 


1.03 


2 1 


4 


2 


0.8 


1.86 


1.03 


2 1 


4 


2 


0.8 


2.5 


3.19 


0.85 


1 2 


2 


0.94 


2.64 


0.97 


1 1 


2 


3 


1.14 


2.11 


1.09 


2 1 


4 


2 


0.8 


2.11 


1.09 


2 1 


4 


2 


0.8 


2.11 


1.09 


2 1 


4 


2 


0.8 


2.11 


1.09 


2 1 


4 


2 


0.8 


2.11 


1.09 


2 1 


4 


2 


0.8 


3 


3.58 


0.88 


1 2 


2 


0.94 


2.99 


1.00 


1 1 


2 


3 


1.14 


2.38 


1.15 


2 1 


4 


2 


0.8 


2.38 


1.15 


2 1 


4 


2 


0.8 


2.38 


1.15 


2 1 


4 


2 


0.8 


2.38 


1.15 


2 1 


4 


2 


0.8 


2.38 


1.15 


2 1 


4 


2 


0.8 


3.5 


3.99 


0.89 


1 2 


2 


0.94 


3.37 


1.03 


1 1 


2 


3 


1.14 


2.67 


1.21 


2 1 


4 


2 


0.8 


2.67 


1.21 


2 1 


4 


2 


0.8 


2.67 


1.21 


2 1 


4 


2 


0.8 


2.67 


1.21 


2 1 


4 


2 


0.8 


2.67 


1.21 


2 1 


4 


2 


0.8 


4 


4.42 


0.91 


1 2 


2 


0.94 


3.76 


1.06 


1 1 


2 


3 


1.14 


2.98 


1.27 


2 1 


4 


2 


0.8 


2.98 


1.27 


2 1 


4 


2 


0.8 


2.98 


1.27 


2 1 


4 


2 


0.8 


2.98 


1.27 


2 1 


4 


2 


0.8 


2.98 


1.27 


2 1 


4 


2 


0.8 


4.5 


4.87 


0.92 


1 2 


2 


0.94 


4.18 


1.08 


1 1 


2 


3 


1.14 


3.31 


1.32 


2 1 


4 


2 


0.8 


3.31 


1.32 


2 1 


4 


2 


0.8 


3.31 


1.32 


2 1 


4 


2 


0.8 


3.31 


1.32 


2 1 


4 


2 


0.8 


3.31 


1.32 


2 1 


4 


2 


0.8 


5 


5.32 


0.93 


1 2 


2 


0.94 


4.61 


1.09 


1 1 


2 


3 


1.14 


3.66 


1.36 


2 1 


4 


2 


0.8 


3.66 


1.36 


2 1 


4 


2 


0.8 


3.64 


1.37 


2 1 


4 


3 


0.98 


3.64 


1.37 


2 1 


4 


3 


0.98 


3.64 


1.37 


2 1 


4 


3 


0.98 


5.5 


5.80 


0.93 


1 2 


2 


0.94 


5.06 


1.11 


1 1 


2 


3 


1.14 


4.02 


1.41 


2 1 


4 


2 


0.8 


4.02 


1.41 


2 1 


4 


2 


0.8 


3.96 


1.43 


2 1 


4 


3 


0.98 


3.96 


1.43 


2 1 


4 


3 


0.98 


3.96 


1.43 


2 1 


4 


3 


0.98 


6 


6.28 


0.94 


1 2 


2 


0.94 


5.52 


1.12 


1 1 


2 


3 


1.14 


4.41 


1.44 


2 1 


4 


2 


0.8 


4.41 


1.44 


2 1 


4 


2 


0.8 


4.29 


1.48 


2 1 


4 


3 


0.98 


4.29 


1.48 


2 1 


4 


3 


0.98 


4.29 


1.48 


2 1 


4 


3 


0.98 


6.5 


6.76 


0.94 


1 2 


2 


0.94 


5.98 


1.13 


2 1 


6 




0.65 


4.81 


1.48 


2 1 


4 


2 


0.8 


4.81 


1.48 


2 1 


4 


2 


0.8 


4.63 


1.54 


2 1 


4 


3 


0.98 


4.63 


1.54 


2 1 


4 


3 


0.98 


4.63 


1.54 


2 1 


4 


3 


0.98 


7 


7.25 


0.94 


1 2 


2 


0.94 


6.35 


1.16 


2 1 


7 




0.71 


5.23 


1.50 


2 1 


4 


2 


0.8 


5.23 


1.50 


2 1 


4 


2 


0.8 


4.98 


1.59 


2 1 


4 


3 


0.98 


4.98 


1.59 


2 1 


4 


3 


0.98 


4.98 


1.59 


2 1 


4 


3 


0.98 


7.5 


7.75 


0.94 


1 2 


2 


0.94 


6.70 


1.20 


2 1 


7 




0.71 


5.66 


1.53 


2 1 


4 


2 


0.8 


5.59 


1.55 


2 1 


5 


2 


0.91 


5.35 


1.64 


2 1 


4 


3 


0.98 


5.35 


1.64 


2 1 


4 


3 


0.98 


5.35 


1.64 


2 1 


4 


3 


0.98 


8 


8.25 


0.94 


1 2 


2 


0.94 


7.07 


1.24 


2 1 


8 




0.76 


6.11 


1.55 


2 1 


4 


2 


0.8 


5.96 


1.60 


2 1 


5 


2 


0.91 


5.72 


1.69 


2 1 


4 


3 


0.98 


5.72 


1.69 


2 1 


4 


3 


0.98 


5.72 


1.69 


2 1 


4 


3 


0.98 


8.5 


8.75 


0.95 


1 2 


2 


0.94 


7.43 


1.28 


2 1 


8 




0.76 


6.57 


1.56 


2 1 


4 


2 


0.8 


6.35 


1.64 


2 1 


5 


2 


0.91 


6.11 


1.74 


2 1 


4 


3 


0.98 


6.11 


1.74 


2 1 


4 


3 


0.98 


6.11 


1.74 


2 1 


4 


3 


0.98 


9 


9.25 


0.95 


1 2 


2 


0.94 


7.80 


1.32 


2 1 


8 




0.76 


7.03 


1.57 


2 1 


4 


2 


0.8 


6.75 


1.68 


2 1 


5 


2 


0.91 


6.50 


1.78 


2 1 


4 


3 


0.98 


6.50 


1.78 


2 1 


4 


3 


0.98 


6.50 


1.78 


2 1 


4 


3 


0.98 


9.5 


9.75 


0.95 


1 2 


2 


0.94 


8.19 


1.35 


2 1 


8 




0.76 


7.51 


1.58 


2 1 


4 


2 


0.8 


7.16 


1.71 


2 1 


5 


2 


0.91 


6.91 


1.82 


2 1 


4 


3 


0.98 


6.89 


1.82 




4 


6 


2.61 


6.89 


1.82 


1 1 


4 


6 


2.61 


10 


10.25 


0.95 


1 2 


2 


0.94 


8.59 


1.38 


2 1 


8 




0.76 


7.99 


1.59 


2 1 


4 


2 


0.8 


7.54 


1.76 


2 1 


7 


2 


1.1 


7.33 


1.85 


2 1 


4 


3 


0.98 


7.25 


1.88 




4 


6 


2.61 


7.24 


1.89 


1 1 


4 


7 


2,88 


10.5 


10.75 


0.95 


1 2 


2 


0.94 


9.00 


1.41 


2 1 


8 




0.76 


8.49 


1.59 


2 1 


4 


2 


0.8 


7.90 


1.82 


2 1 


7 


2 


1.1 


7.76 


1.88 


2 1 


4 


3 


0.98 


7.61 


1.94 




4 


6 


2.61 


7.60 


1.95 


1 1 


4 


7 


2.88 


11 


11.25 


0.95 


1 2 


2 


0.94 


9.43 


1.44 


2 1 


8 




0.76 


8.98 


1.59 


2 1 


4 


2 


0.8 


8.26 


1.88 


2 1 


8 


2 


1.18 


8.16 


1.93 


1 1 


4 


5 


2.31 


7.99 


2.00 


1 1 


4 


6 


2.61 


7.96 


2.02 


1 1 


4 


7 


2.88 


11.5 


11.75 


0.95 


1 2 


2 


0.94 


9.87 


1.46 


2 1 


8 




0.76 


9.45 


1.61 


1 1 


4 


3 


1.86 


8.62 


1.94 


2 1 


8 


2 


1.18 


8.55 


1.97 


1 1 


4 


5 


2.31 


8.36 


2.06 




4 


6 


2.61 


8.32 


2.08 


1 1 


4 


7 


2.88 


12 


12.25 


0.95 


1 2 


2 


0.94 


10.32 


1.47 


2 1 


8 




0.76 


9.90 


1.62 


1 1 


4 


3 


1.66 


8.99 


2.00 


2 1 


8 




1.18 


8.95 


2.02 


1 1 


4 


5 


2.31 


8.75 


2.12 




4 


6 


2.61 


8.70 


2.14 


1 1 


4 


7 


2.88 


12.5 


12.75 


0.95 


1 2 


2 


0.94 


10.77 


1.49 


2 1 


8 




0.76 


10.37 


1.63 


1 1 


4 


3 


1.66 


9.37 


2.06 


2 1 


8 


2 


1.18 


9.36 


2.06 


1 1 


4 


5 


2.31 


9.14 


2.17 




4 


6 


2.61 


9.07 


2.20 


1 1 


4 


7 


2,88 


13 


13.25 


0.95 


1 2 


2 


0.94 


11.24 


1.50 


2 1 


8 




0.76 


10.82 


1.65 


1 1 


8 


2 


2.02 


9.76 


2.11 


2 1 


8 


2 


1.18 


9.76 


2.11 


2 1 


8 


2 


1.18 


9.52 


2.23 


3 1 


9 


2 


0.84 


9.46 


2.26 


1 1 


4 


7 


2,88 


13.5 


13.75 


0.95 


1 2 


2 


0.94 


11.72 


1.51 


2 1 


8 




0.76 


11.19 


1.70 


1 1 


8 


2 


2.02 


10.17 


2.16 


2 1 


8 


2 


1.18 


10.17 


2.16 


2 1 


8 


2 


1.18 


9.91 


2.28 


3 1 


9 


2 


0.84 


9.84 


2.32 


1 1 


4 


7 


2,88 


14 


14.25 


0.95 


1 2 


2 


0.94 


12.20 


1.51 


2 1 


8 




0.76 


11.57 


1.75 


1 1 


8 


2 


2.02 


10.58 


2.20 


2 1 


8 


2 


1.18 


10.58 


2.20 


2 1 


8 


2 


1.18 


10.28 


2.36 


2 1 


8 


3 


1.52 


10.23 


2.38 


1 1 


4 


7 


2.88 


14.5 


14.75 


0.95 


1 2 


2 


0.94 


12.68 


1.52 


2 1 


8 




0.76 


11.97 


1.79 


1 1 


8 


2 


2.02 


11.01 


2.24 


2 1 


8 


2 


1.18 


10.99 


2.25 


2 1 


9 


2 


1.25 


10.65 


2.43 


2 1 


8 


3 


1.52 


10.63 


2.44 


1 1 


4 


7 


2.88 


15 


15.25 


0.95 


1 2 


2 


0.94 


13.17 


1.52 


2 1 


8 




0.76 


12.38 


1.83 


1 1 


8 


2 


£.02 


11.45 


2.27 


2 1 


8 


2 


1.18 


11.38 


2.30 


2 1 


10 


2 


1.35 


11.02 


2.50 


2 1 


8 


3 


1.52 


11.02 


2.50 


2 1 


8 


3 


1.52 



TABLE IV 



Table of the P-CPM schemes with best Cp_cpm and raised-cosine frequency pulse at 99% bandwidth. The complexity range is 8 
TO 512. The highlighted entries correspond to the selected CPM SCHEMES for target CAPACITIES OF 1, 1.5 AND 2 bits/s/Hz. 




Fig. 7. The optimized CPM encoder. The block labelled Ep (resp. Sm) is 
a modulo-P (resp. modulo- Af) adder. The block labelled G _1 is the inverse 
of the Gray function defined in [18, ch.20]. 



schemes with best Cp- 
and Tab. |IV] 



CPM- 



Results are shown in Tab. 



Ill 



A comparison of Tab. [EITJwith Tab. [I] and Tab. [IV] with Tab. [TT] 
highlights that, in some cases, the CPM schemes with best 
Ccpm feature the best Cp_cpm. as is the case, e.g., for REC 
pulses with complexity 64 and capacity 1.5 bits/s/Hz and for 
RC pulses with complexity 64 and capacity 1.5 bits/s/Hz and 
2 bits/s/Hz. In some other cases, the selected schemes do not 
match. 

We observed that, in many binary CPM schemes, the 
pragmatic capacity resulting from the CPE optimization ap- 
proaches the CPM capacity closer than for non binary 
schemes. As a result, even when featuring a worse CPM ca- 
pacity, these schemes may exhibit a larger pragmatic capacity. 
This is the case, e.g., of REC pulses with complexity 64 and 
capacity 1 bits/s/Hz, where the scheme with best Ccpm is 
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quaternary and the scheme with best Cp_cpm is binary. 

The optimal mapping described above results in a CPM 
scheme with finite d^i, the minimum Euclidean distance at 
the CPM output corresponding to an input Hamming distance 
1. To prove this, we use the polynomial notation with dummy 
variable D for symbol sequences and we observe from Fig. [7] 
that 



fli+ip) = Doi[p), 



1. 



1 



and 



a L+1 (D) = a L (D) 



D 

1 -D 



(V.19) 



(V.20) 



If Q = 1, from (V20) , ( |VT9) and ( fVT8| , we obtain 

a 1 ( J D) = (i(£»)(l-D)) M . 
Therefore ai(-D) can be obtained by feedforward encoding 



of i{D). From (V.19 1 it is straightforward to deduct that all 

a,i(D), i = 1, 



, L have the same property. 
As for symbol ql+i{D), from ( V.20[ ) we obtain 



a L +i(D) = 



D L ai (D) 
1 - D 



D L (i(D)(l-D)) M 
1 — D 



(Vil) 



If M = zP, z 6 N, this clearly proves that symbol cll + i(D) 
can be obtained by feedforward encoding i(D). When P = 
wM, w € N, it is possible to rewrite ( V.21 1 as 



a L+ i(L>) 



D L (i(D)(l-D)) 



M 



M 



1 - D 
D L (i(D)(l-D)) M 
M(l - Z>) 

(^^(^))m 

(*(£>)(!-£>)) 



M 



P/M 



M 



M 



M(l - D) 



(V.22) 



P/M 



When the last term of (V.22 1 is zero, et£ + i(Z3) can be obtained 



by feedforward encoding from i(D), therefore the encoder 
is feedforward. For difference sequences i(D) with length 1 
and symbols in {— 1,+1} this condition is met. Moreover, 
due to Gray mapping, the corresponding pairs of binary input 
sequences have unitary Hamming distance, hence the resulting 
encoder has finite dei- 

It is well known that in serially concatenated schemes, in 
order to take advantage of the interleaver gain [19], the inner 
constituent encoder or modulator must have infinite dsi [17]. 
For this reason, when used in SC-CPM schemes, the optimal 
mapping leads to poor performance and hence a recursive CPE 
must be employed. We have thus used the mapping induced 
by the Rimoldi decomposition, which has this property. 

VI. Implementation of the optimized coded CPM 

SCHEMES 

Using the procedures described in the previous sections, 
SC-CPM and P-CPM schemes have been designed with target 
capacities of Ct = 1 bits/s/Hz, Ct = 1.5 bits/s/Hz and Ct = 
2 bits/s/Hz. 



A. Selection of the modulation schemes 

Fig. [8] and Fig. [9] show the Ei,/Nq values plotted against the 
CPM complexity for different target capacities. Continuous 
curves refer to CPM capacity and dashed curves refer to 
pragmatic capacity. 

For REC schemes (Fig. [HJ we observe that increasing the 
complexity beyond 64 yields a negligible gain both in CPM 
capacity and in pragmatic capacity. Hence, CPM schemes with 
REC frequency pulse and maximum complexity 64 have been 
considered. 

For RC schemes (Fig. |5J, increasing the complexity beyond 
64 yields a significant gain both in CPM capacity and in 
pragmatic capacity at 1.5 bits/s/Hz and 2 bits/s/Hz, hence CPM 
schemes with higher complexities should be chosen. However, 
in SC-CPM receivers, the increased complexity would be 
enhanced by the iterative process, resulting in an impractical 
solution. 



- 2 bits/s/Hz (P-CPM) 
-2 bits/s/Hz (CPM) 

- 1.5 bits/s/Hz (P-CPM) 
-1.5 bits/s/Hz (CPM) 

- 1 bits/s/Hz (P-CPM) 
- 1 bits/s/Hz (CPM) 




Fig. 



CPM complexity 

Energy efficiency of REC CPM schemes. 



■ 2 bits/s/Hz (P-CPM) 
-2 bits/s/Hz (CPM) 
-1.5 bits/s/Hz (P-CPM) 

- 1 .5 bits/s/Hz (CPM) 

- 1 bits/s/Hz (P-CPM) 
1 bits/s/Hz (CPM) 




CPM complexity 

Fig. 9. Energy efficiency of RC CPM schemes. 



B. Serially-concatenated schemes 

A serially concatenated CPM scheme consists of the cascade 
of an outer convolutional encoder connected to the CPM 
modulator through an interleaver. 

As outer encoder, a 4-state, rate 1/2 systematic recursive 



convolutional encoder has been chosen. It is shown in Fig. 10 
The choice has been dictated by a trade-off between perfor- 
mance and complexity. The encoder output can be punctured 
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Fig. 10. Convolutional encoder used in SC-CPM schemes. 



Fig. 11. Scheme of the CPM modulator with m = 2, h = 1/5, L = 2. 



to achieve rates up to 3/4 while keeping a free distance larger 
than 2. 

In order to obtain rates larger than 1/2, the output of 
the convolutional encoder is punctured according to a rate- 
matching algorithm that selects all systematic bits and some 
coded bits. The puncturing rate is 



Rp 



No 
Nj 



> 1 



where N and TV/ are two integers (Nj < N < 2TV/). The 
puncturing algorithm selects Nj coded bits from each block 
of Nq with the following constraint: all systematic bits are 
selected; some coded bits are selected to achieve the desired 
rate; the other bits are punctured. The resulting outer code rate 



Rcc — 



INo 
2 TV/ 



The interleaver that connects the punctured outer encoder 
to the inner CPM modulator is a spread interleaver. 

The design of a SC-CPM scheme consists in the choice 
of the CPM scheme and of the outer encoder based on the 
desired capacity and on design constraints (e.g., the maximum 
affordable complexity of the CPM scheme). 

Example 6.1: We describe the design of a SC-CPM scheme 
with target capacity Ct = 1.5 bits/s/Hz, REC frequency pulse 
and maximum CPM complexity of 64 edges per CPM input 
bit. As for SC-CPM schemes, the parameters that maximize 
Ccpm must be chosen, from Tab. [I] Moreover, since the 
maximum CPM complexity is 64, we must restrict the search 
to the column set labelled 64. Next, in the column labelled C, 
we look for the value closest to the target capacity. We find that 
the capacity value closest to Ct is 1.51 and the corresponding 
CPM parameters are m = 2, h = 1/5, L = 2. Such scheme 
exhibits a symbol rate R s = 1.18 and a corresponding bit 
rate Rb = mR s = 2.36. The scheme of the CPM modulator 



corresponding to these parameters is shown in Fig. 1 1 

In order to obtain a target capacity Ct = 1.5 bits/s/Hz, the 
outer code rate must be set to 



therefore 



Rcc 



No 
TV/ 



2Nj 



2C 7 



Ct 
Ri, 



0.636 



Rh 



1.2712 



75 
59' 



The final step consists in choosing two integers for No and 
Nj in order to approximate with sufficient precision this ratio. 

□ 

The resulting design parameters for all considered target 
capacities are shown in Tab. [V] 

C. Pragmatic schemes 

A pragmatic CPM scheme consists of the cascade of a 
SCCC encoder and a CPM modulator. 

The chosen encoder is the SCCC encoder described in [12]. 
The constituent convolutional encoders used in such scheme 
are the 4-state convolutional encoders shown in Fig. 



10 The 



outer convolutional encoder is punctured to a rate 2/3. 

The inner encoder is punctured using a carefully designed 
rate-matching algorithm that selects all the systematic bits and 
some coded bits in order to achieve the desired overall rate. 

The interleaver that connects the punctured outer encoder 
to the inner encoder is a spread interleaver. 

Example 6.2: We show how to design a P-CPM scheme 
with target capacity of 2 bits/s/Hz, RC frequency pulse and 
maximum CPM complexity of 64 edges per CPM input bit. For 
P-CPM schemes, the inner CPM modulation parameters that 
maximize Cp_cpm must be chosen. Therefore, Tab. IIVI must 



be considered. Moreover, since the maximum CPM complexity 
is 64, we must restrict the search to the column set labelled 64. 
Next, in column labelled C we look for the value closest to the 
target capacity. We find that the capacity value closest to Ct is 
2.00 and the corresponding CPM parameters are to = 2, h = 
1/8, L = 2. Such scheme exhibits a symbol rate R s = 1.18 
and a corresponding bit rate Ri, — mR s — 2.36. The scheme 
of the CPM modulator corresponding to these parameters is 



shown in Fig. 12 




Fig. 12. Scheme of the CPM modulator with m = 2, h = 1/8, L = 2. 

In order to obtain a target capacity Ct = 2 bits/s/Hz, the 
outer code rate must be 
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Ct 




SC-CPM 






P-CPM 






2.5 




A/V-, 


Nj 


m 


p 




^SCCC 


m 








u 


i 
i 


ZJ 


24 


2 




Z 


0.752 


1 


A 






2 


1.5 


75 


59 


2 


5 


2 


0.781 


2 


4 


2 


2 


2 


40 


27 


2 


6 


2 


0.735 


1 


4 


4 






1 


5 


4 


2 


4 


2 


0.625 


2 


4 


2 




u 


1.5 


150 


91 


2 


5 


2 


0.938 


2 


4 


2 


X 

^ 1.5 




2 


100 


59 


2 


8 


2 


0.847 


2 


8 


2 


.5 



TABLE V 

Parameters of the optimized SC-CPM and P-CPM schemes 

YIELDING 1.0, 1.5 AND 2.0 BITS/S/HZ. 



i?sccc = = 0.847. 
Kb 



□ 

Tab. [V] summarizes the parameters of the designed schemes. 

VII. Results 

The performance of the designed coding schemes reported 
in Tab. [V] has been assessed through simulation. For each 
scheme, the CPM capacity and the pragmatic capacity have 



been estimated. The obtained results are shown in Fig. 13 for 



RFC frequency pulses and Fig. 14 for RC frequency pulses 




E b /N„ [dB] 

Fig. 13. CPM capacity and pragmatic capacity of the selected CPM schemes 
with complexity y < 64. The frequency pulse is REC. 



In both the P-CPM and the SC-CPM case, an outer infor- 
mation word length of K = 9840 bits has been considered 
and the code rates have been set according to Tab. [V] The 
number of decoding iterations has been set to 10. 

The information rates Tsc-cpm (resp. -Tp-cpm) of the 
SC-CPM (resp. P-CPM) schemes has been evaluated through 
simulation by computing the mutual information between the 
information bits at the input of the channel encoder and the 



soft outputs of the channel decoder. Fig. 15 Fig. 16 and Fig. 17 
show the obtained information rates for RFC frequency pulse 
and complexity y < 64. We observe that, for a capacity 
of 1 bits/s/Hz, the pragmatic scheme achieves the desired 
information rate at a lower E), /N n than the SC-CPM scheme, 
and at 1 dB from the pragmatic capacity. For capacities of 1 .5 
bits/s/Hz and 2 bits/s/Hz, the SC-CPM scheme achieves the 




E b /N„ [dB] 

Fig. 14. CPM capacity and pragmatic capacity of the selected CPM schemes 
with complexity 3^ < 64. The frequency pulse is RC. 



desired information rate at a lower E^/Nq than the P-CPM 
scheme. 

The best SC-CPM and P-CPM for the chosen spectral 
efficiencies have been simulated over an AWGN channel so 
as to obtain accurate estimates of the bit and frame error 
probabilities at very low values. The results are shown in 
We observe that, at 1 bits/s/Hz, the P-CPM scheme, 



Fig. 18 



while featuring a lower complexity (see Tab. [Vfl, exhibits a 
0.5 dB gain over the SC-CPM scheme at FER = 10~ 4 . At 
higher spectral efficiencies, the SC-CPM scheme exhibits a 
gain of more than 1 dB over the P-CPM scheme. 



Fig. 19 shows similar results for RC frequency pulses. The 
SC-CPM schemes feature a 0.3 dB gain over the P-CPM 
schemes at 1 bits/s/Hz and 2 bits/s/Hz, while at 1.5 bits/s/Hz 
the gain of the SC-CPM scheme is about 1 dB. In both 
frequency pulses and all spectral efficiencies, no error floors 
have been observed. 




E b /N„ [dB] 

Fig. 15. Information rates of P-CPM and SC-CPM coded systems at 1 
bits/s/Hz. y = 64, REC frequency pulse. 



A comparison of REC and RC pragmatic schemes is shown 



in Fig. 20 We observe that the RC scheme performs better 
than REC scheme at 1 bits/s/Hz and 1.5 bits/s/Hz. At 2 
bits/s/Hz, the hierarchy is reversed. At 2 bits/s/Hz, REC 
schemes feature a gain of roughly 0.6 dB over the RC schemes 
at FER = 10~ 4 . 
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E b /N [dB] 

Fig. 16. Information rates of P-CPM and SC-CPM coded systems at 1.5 
bits/s/Hz. y = 64, REC frequency pulse. 
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Fig. 19. Error rates of P-CPM (curves without markers) and SC-CPM (curves 
with markers) coded systems, y = 64, RC frequency pulse. 
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Fig. 17. 
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y = 64, REC frequency pulse. 



Fig. 20. Comparison of REC (curves without markers) and RC (curves with 
Information rates of P-CPM and SC-CPM coded systems at 2 markers) P-CPM schemes with y < 64. 
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Fig. 18. Error rates of P-CPM (curves without markers) and SC-CPM (curves 
with markers) coded systems, y = 64, REC frequency pulse. 



In order to perform a fair comparison, the decoding com- 



plexity must be kept into account. Using definitions (II. 3 i and 



(II.4i, the decoding complexities of the two schemes have 
been computed. For the P-CPM scheme, the complexity of the 
MHOMS [12] binary decoder has been computed as follows: 



ccc 



= 2N lt N so + -N s 



(VII. 1) 



where the factor 3/2 derives from the outer code rate of the 
MHOMS scheme, which is 2/3. Results are shown in Tab. [VT| 
we observe that the complexity of the SC-CPM scheme is 
always larger than that of the P-CPM scheme. We define 
the complexity ratio Ry = IVsc-cpmA'p-cpm: its value 
is larger than 2 and exceeds 3 for a 1 bits/s/Hz capacity 
with REC frequency pulse. Higher coding gains could be 
achieved increasing the P-CPM complexity, thus reducing the 
performance gap with respect to SC-CPM schemes. 





Capacity 


P - CPA 

ysccc ycpM 


y p-cpm 


SC - CPM 

y SC-CPM 


Ratio 
Ry 


REC 


1 

1.5 

2 


200 21 
200 41 
200 87 


221 
241 
287 


694 
709 
728 


3.14 
2.94 
2.54 


U 


1 

1.5 

2 


200 51 
200 34 
200 76 


251 
234 
276 


592 
565 
733 


2.36 
2.41 
2.66 



TABLE VI 

Complexity of P-CPM and SC-CPM schemes, y - 



VIII. Conclusions 

The pragmatic approach to coded continuous-phase modula- 
tion (CPM) has been proposed as a capacity-approaching low- 
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complexity alternative to the serially-concatenated CPM (SC- 
CPM) coding scheme. After performing a selection of the best 
spectrally-efficient CPM modulations to be embedded into SC- 
CPM schemes, the pragmatic capacity of CPM modulations 
has been evaluated and optimized through a careful design 
of the mapping between input bits and CPM waveforms. 
The so obtained schemes have been cascaded with an outer 
serially-concatenated convolutional code to form a pragmatic 
coded-modulation system. The resulting schemes have been 
shown to exhibit performance close to the CPM capacity 
without requiring iterations between the outer decoder and the 
CPM demodulator. As a result, the receiver exhibits reduced 
complexity and increased flexibility due to the separation of 
the demodulation and decoding functions. 
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Appendix I 
Proof of Theorem [5J] 

First, we revise some properties of the CPM trellis, for a 
binary CPM scheme. 

An error event generated by a given difference sequence 
b = (bx, ■ ■ ■ , b/\) £ {— 1, 0, 1} A , is a pair of trellis paths with 
length A + L — 1 trellis steps, which diverge at, say, time 
and merge back together at time A + L — 1. For a given value 
of b, there exist several different error events, corresponding 
to all possible choices of the variables: 



Po,a- 



L+2, 



a „W 

,a ,a 1 ,. 



and 



Po, o,-l+2, ... ,ao,a x 



(i) 



(2) 

, O a , QA+1, 



•j a A +L-i (LI) 



, OA+L-l (L2) 



where f3 ,ci- L ^ 



, ap determine the starting state, 



and 



,(2) = 




(1.3) 



(1.4) 



and oa+i, . . . , ga+l-i determine the final state of the error 
event. If two trellis edges belong to the same trellis section of 
the error event, we call them an edge pair of the error event. 

Now, we build the graph G(b) = (V, £) as described 
in Sect. 



V-B 



two vertices are adjacent if and only if the 
corresponding trellis edges do not have the same starting state 
and constitute an edge pair of an error event generated by b. 
The graph G(b) has Af(b) connected components, denoted 
Co, . . . ,C/v"(b)-i> where C{ = (Vi,£i), whose properties will 
be investigated hereafter. 

For an integer amount 7, the ^-rotation transforms the trellis 
edge (a, (3) into the trellis edge (a, ((3 + y)p). The image of 
a set of edges through the 7-rotation is the set of images of 
the edges belonging to the set through the 7-rotation. The 
following proposition holds: 

Proposition 1.1: For any 7, the 7-rotation transforms the 
set Vi, corresponding to a component C, of G(b), into a set 
Vj, corresponding to another component Cj (even, possibly, 
the same). 

Proof: Two trellis edges (a, f3) and (a' , f3') are in the 
same component Cj if and only if in the graph there is a length- 
re path, 



(a (K) ,/3 (K) ) 
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with (a«,/?W) = (a,/3) and (a^,^) = (a',/3')- 
This means that, for i = 1, ,..,n — 1, (cr 8 -*,/?") and 
have different starting states and form an edge 
pair of an error event generated by b. 

Now, notice that, given an error event generated by b, by 



changing /3q in (1.1 1 and ( 1.2 1 into (/?o+7) P> we obtain another 



error event generated by b, whose edge pairs are in one-to-one 
correspondence through 7-rotation with the edge pairs of the 
original error event. For this reason, in the graph there is the 
path 



(a« /3«)-(a< 2 >,/f>) 



(a (K) ,M K) ) 



where ffl ^ (0W+ 7 ) P . 

Thus, by definition, (a, (f3+-f)p) and (a', (fi'+'fjp) belong 
to the same component Cj. For the arbitrariness of (a,/?) and 
we have that the image of Vj through 7-rotation is 
contained in Vj. 

Finally, notice that, with the same argument, the image of 
Vj through (P — 7) -rotation is contained in Vi- From this fact, 
and the injectivity of rotations, we deduce that the image of 
Vi through 7-rotation is equal to Vj. ■ 

From (I.li and ( |I.2| i, we can deduce a set of rules for a pair 
of trellis edges to be in the same component, in the following 
way: 

• We fix a trellis section inside the error event. The error 
event has length A + L — 1 trellis sections, but the first 
is not to be considered, thus we have A + L — 2 possible 
choices. 

• We consider the constraints imposed on an edge pair of 
the error event at that trellis section. 

We derive two types of rules, according to the trellis section 
we are considering. 

Type-I rules: They correspond to the last L trellis sections 
of the error event. They give the following edge pairs, for 
j=A-L+l,...,A: 



A+L-l 



and 



(»f. (0 



U (*) I W 

where ay = [ay, 
and 



A+L-l 



,a% L _^af = aj if j < 1 or j > A, 



- Qwh (af 
Qw H (af) 



(1.5) 



(1.6) 



/0A+L-1 = A) + Q 



3 
A 

;=-L+2 



(i) 



iA (1) 



E 



A 



(2) 



The above relation holds because 
mod P. 

Type-II rules: They correspond to the A — 2 trellis sections 
of the error event from the second to the (A — l)-th one. They 
give the following edge pairs, for j = —L + 3, . . . , A — L: 



(i) 0(1) 



j+L-l 



and 



a 



(2) fl (2) 



j+L-l 



(1.7) 
(1.8) 



where a® = (a^\ . . . , o^ i _ 1 ), aj 1 - 1 = aj if j < 1 or j > A, 
and 



(i) 



Notice that, if A = 2, there are no type-II rules. 

The following lemma will be the key to prove Theorem |5T 



Lemma 1.2: Consider the sets of trellis edges: 

V t = {(a, (i - Qw H {a)) P ) : a € {0, 1} L } 



(1-9) 



for i = 0, . . . , P — 1. Each of these sets is entirely contained 
in one connected component of the graph 5(b). 

Proof: We prove by induction that the set of type-I rules 
implies the lemma, starting from the last rule (i.e., j = A in 



(1.5 1 and (1.6 1) and going backwards. 

Suppose that the last d rules (d < L) imply that the set of 
edges of the form: 



V / 3 A+£ _ 1 , d (x) = {([d,x] ! (^ AH 



L-l 



Q^([d,x])) p ),de{0,i} d } 



for a fixed length-(L — d) vector x = (xi, . . . ,xl-<i), is 
entirely contained in one connected component. This is true 
for d = 1, since b& ^ 0. 

Now, again because &a 7^ 0, the last-but- (d + 1) rule 
pairs edges like d'xi . . . a^-d and d"xi . . . XL-d, for suitable 
length-G? vectors d' and d", where x is the complement of x. 
This implies that 



Vp A+L _ u d+i{x) = V/3 A+£ _ 1 , d ([0,x]) U Vft, 



i([l,x]), 



for a fixed length-(L — d — 1) vector x, is entirely contained 
in one connected component. 

By induction up to d = L, we have proved that, for all 
possible values of /?a+l-i, = Vp A+t ,_ u L is entirely 

contained in one connected component. ■ 
Proof: [Of Theorem |5T] If A(b) = 2, then there are 
only type-I rules, and the components of 5(b) satisfy, for 
i = 0,...,P-l: 

H = Vi (1. 10) 



and the theorem is proved. 

If A(b) > 2, instead, there are also type-II rules. Consider, 
in particular, the type-II rule with j = A — L in (1.7 1 and (1.8 1. 
It gives the edge pair: 



3 {1) 

, "A-L'PA-1 



and 



(2) 



(2) 
A-l 



(1.11) 



(1.12) 



and to 



But the above transitions belong to V( | g A+£ _ 1 _g) I 
V/3 A+L _ 1 . Then, these two sets are merged by the type-II rule 
with j = A — L. Consider the component C\ of 5(b) that 
contains these sets and let Vi be its vertex set. Prop. |1.1| 
implies that a Q-rotation transforms V\ into itself. But then, 
C\ also contains V( ( 3 a+l _ 1 _;q) p , for every integer /. Since Q 
and P are 3, we deduce that C\ contains all Vj's and thus is 
the only component of 5(b). This concludes the proof. ■ 
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Appendix II 
Proof of TheoremL 

To prove Theorem 5.2 we first notice that Prop. 1 1 . 1 1 holds 
for a 1 CPM scheme also, since in the proof, the hypothesis 
of binary input symbols is not used. 

Also, there are two types of rules, which are the extension to 
the 1 case of the two types introduced in the previous sections. 

Type-I rules: They correspond to the last L trellis sections 
of the error event. They give the following edge pairs, for 
j=A-L 

L-l 



and 




(II. 1) 



(H.2) 



where ay = {a^\ . . . , a^ i _ 1 ), — a 3 if j < 1 or j > A, 
and 

\ i=-L+2 / p 

The above relation holds because 2~2iLi a i ~ 2~2iLi 
mod P. 

Type-II rules: They correspond to the A — 2 trellis sections 
of the error event from the second to the (A — l)-th one. They 
give the following edge pairs, for j = —L + 3, . . . , A — L: 



and 



where = (a^\ . . . , a^] rL _ l ), = a 3 if j < 1 or j > A, 
and 

\ k=-L+2 ) p 

Notice that, if A = 2, there are no type-II rules. 

The following lemma is the extension of Lemma 1.2 to the 
1 case. 

Lemma 2.1: Let b be a difference sequence of length A 
with 6a = ±1- Consider the sets: 

5 = 1^, (i-Qj2 a *J ) :ae{0,...,M-l} L j, 

(115) 

for i = 0, . . . , P — 1. Each of these sets is entirely contained 
in one connected component of the graph £7(b). 

Proof: By induction on the set of type-I rules. It is the 
straightforward extension of Lem ma |1.2| ■ 
Finally, the proof of Theorem |5.2| follows the proof of 
Theorem 5.1 by noticing that if A(b) = 2, then there are 
only type-I rules, while if A(b) > 2 there are also type-II 
rules, that cause the graph Q(b) to be connected. 



